Properties of a coupled two species atom-heteronuclear molecule condensate 
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We study the coherent association of a two species atomic condensate into a condensate of het- 
eronuclear diatomic molecules, using both a semiclassical treatment and a quantum mechanical 
approach. The differences and connections between the two approaches are examined. We show 
that, in this coupled nonlinear atom-molecule system, the population difference between the two 
atomic species play significant roles in the ground state stability properties as well as in coherent 
population oscillation dynamics. 

PACS numbers: 03.75.-b, 05.30.Jp 



I. INTRODUCTION 

After the experimental realization of the trapped 
atomic Bose- Einstein condensates (BECs), achieving 
molecular BEC has been regarded as another milestone 
in the field of ultracold atomic physics. As molecules 
are inherently much more complex in energy spectrum 
than their constitutes-atoms, direct laser cooling meth- 
ods popular with atoms are ineffective with molecules. 
Much recent activities, both in experiments Q], 3, |3, 
1 [ESS and in theory & M, M, 03, lii 

snaiii os m mi, hi in, m, have been focused 

primarily on converting ultracold atoms into ultracold 
molecules by means of magneto- (Feshbach resonance) 
or photo-association, in which two atoms are combined 
into a diatomic molecule mediated by either a mag- 
netic field or an optical field. Both ultracold degenerate 
bosonic and fermionic atoms have been successfully con- 
verted into molecules. Considerable theoretical efforts 
have been devoted to improving the conversion efficiency 
[1, M, E3, 11 , H3 and understandin g th e molecular as- 
sociation pi, Gj, E M, EE M, Erp, p, rjp. M, HI 

as well as the dissociation dynamics [25l [2g, [27| of the 
atom-molecule coupling model. 

It needs to be emphasized that most of the aforemen- 
tioned studies, with the notable exception of Refs. [j3] 
and [HI, concern homonuclear molecules. The interest of 
this paper is, however, the heteronuclear molecules in the 
coupled atom-molecule systems with two different atomic 
species. As a natural progression, quantum degener- 
ate heteronuclear molecules are expected to be the next 
challenge to the atomic physics community, because het- 
eronuclear molecules possess intriguing properties that 
will open up many new avenues of research. For example, 
unlike their homonuclear counterpart which are always 
bosonic, heteronuclear diatomic molecules can be cither 
bosons or fermions, hence quantum statistics will play 



important roles in such systems [1 21 ] . Furthermore, large 
electric dipole moment can be induced in heteronuclear 
molecules with the prospect of creating dipolar superfluid 
[28| and with potential applications in quantum comput- 
ing [29( | , quantum simulation (3fj| and test of fundamental 
symmetry [3ll ] . For these reason, heteronuclear molecules 
have recently received much theoretical and experimental 
attention. Already, Feshbach resonances have been ob- 
served in various quantum degenerate Bose-Fermi atomic 
mixtures [H, H, S| , and heteronuclear molecules from both 
Bose-Fermi and Bose-Bose mixtures have beenproduced 
through the photoassociation technique [4j, |5|, |6( . 

In this paper, we consider, within a three-mode model, 
a system of bosonic diatomic heternuclear molecules cou- 
pled to its constituent atoms, both types of which are also 
assumed to be bosonic. Besides the collisional strengths 
and the detuning (bare energy difference between the 
molecular and the atomic modes), due to the presence 
of two types of atoms, we have a new "control knob" 
— the population imbalance between the two species — 
which we shall pay special attention to. We note in pass- 
ing that recent experiments on two-component degener- 
ate Fermi gases with population imbalance [32l . 1 3 3l ] have 
generated great excitement due to its rich phase diagram 
with various exotic quantum phases in which the popu- 
lation imbalance plays a critical role. We will study our 
system using both a mean-field semiclassical and a full 
quantum mechanical method. The differences as well as 
the connections between the two approaches will be ex- 
amined. 

The paper is organized as follows. In Sec. II we present 
our model in both full quantum and the mean-field ver- 
sion. In Sec. Ill we study the ground state properties and 
their relevance in creating the molecules from the atoms 
by adiabatically sweeping the detuning. The population 
dynamics is presented in Sec. IV and finally we con- 
clude in Sec. V. Our work differs from Refs. [12] and [l3j] 
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in the following ways: Ref. p3| focuses on the quantum 
statistical properties of the molecules and does not con- 
sider the effect of population imbalance; while Ref. [HI 
uses a very different quantum approach (Bethe ansatz) 
from ours and does not pay much attention to the atom- 
molecule conversion process. 



II. QUANTUM MODEL AND MEAN-FIELD 
APPROXIMATION 



We adopt a simple three-mode model in which we de- 
scribe our atom-molecule system with two atomic modes 
(1 and 2) and one molecular mode (m). The basic as- 
sumption here is that the spatial wave functions for these 
modes are fixed so that we can associate each mode 
with an annihilation operator hi of a particle in mode 
i{= 1,2 and m). Similar models have been extensively 
used in the studies of condensates in double-well po- 
tentials 1341 . |47L [481 ] ; coupled atom-molecule condensates 
[13, ED, EI ISH lH, EE E3, E d, as well as spinor 
condensates [35j . 

Within the three-mode approximation, the second 
quantized Hamiltonian reads 



H = S (V m a m +g {h} m hih 2 + h.c.) + 



^2xijala]ajai , (1) 



where the detuning 8 represents the energy difference 
between the molecular and atomic levels which can be 
tuned by external field, g is the atom-molecule coupling 
strength and Xij = Xji is the s-wave collisional strength 
between modes i and j. Without the collisional terms 
our model will reduce to the trilincar Hamiltonian de- 
scribing the nondegenerate parametric down-conversion 
in quantum optics [36l [37| • 

There are two obvious constants of motion from Hamil- 
tonian (TT]): 



N = a\ai + h\h 2 + Sh^h, 



D 



a\ai - a) 2 h 2 , (2) 



which account for the total particle number and the num- 
ber difference between the two atomic species, respec- 
tively. Taking advantage of the constants of motion, the 
Hamiltonian (fT]) can be simplified as 



H 




(a\ai + h\h 2 



(3) 



where we have introduced two dimcnsionlcss quantities 

A = N (Xll + X22 + X mm + 2 Xl2 - 2 Xml ~ 2 Xm2) /G , 

A = [8-(D-l) Xll + (D + l) X22 + (N-l) Xmm 
-(N-D) Xml -(N + D) Xm2 ]/G, 

with G = g\/2N as the rescaled atom-molecule coupling 
strength. In writing ([3]), we have neglected the constant 
terms proportional to D and N. 



To complement the quantum study, we develop a semi- 
classical description of our system by following the usual 
mean-field approach, which has proven to be a powerful 
tool for the study of Bose-Einstein condensates. As a first 
step, we apply the Heisenberg equation to arrive at the 
operator equation for hi and then replace hi with the cor- 
responding c- number aj. Next, we change the equation 
for a$ into the ones for Ni and ip i through the transfor- 
mation di = \fWie l ' Pi , where AT, and ipi represent the 
number and phase of the bosonic field for the particles 
in species i, respectively. Finally, we take advantage of 
the existence of the two conserved quantities N and D, 
and simplify our problem into a one described by two 
variables: the normalized population in the two atomic 
modes 

x = (N 1 + N 2 ) /N , 

and the phase difference 

<p = <p x + <p 2 - tp m . 

The equations of motion for x and tp can be easily ob- 
tained as 



— — = — \/ (1 — x) (x 2 — cP) sin tp , 
or 



A- Ax 



d 2 + 2x - 3x 2 

2y/(l- X ) (x 2 -d 2 ) 



(4a) 

cos (p , (4b) 



where r = Gt is the dimensionless time and d — D /N the 
normalized atomic population imbalance. Without loss 
of generality, we will assume a non- negative d € [0, 1]. 

In the language of Hamiltonian mechanics, x and ip 
form a pair of canonically conjugate variables satisfying 
the equations 



dx 



dH 

dip 



dip 



dH 

dx 



with the dimensionless mean-field Hamiltonian H given 
by 



H 



A 



Ax + y/(l-x) (x 2 -d 2 )cosp. (5) 



We note that if d = 0, i.e., when the two atomic modes 
have the same population, Hamiltonian ([5]) would have 
the same form as the corresponding Hamiltonian de- 
scribing homonuclear molecule association from a single 
atomic mode flU [l{| . The quantum mechanical Hamil- 
tonian and its semiclassical counterpart serve as 
the starting point of our study. 



III. STEADY STATES AND RAPID ADIABATIC 
PASSAGE 

Semiclassically, the fixed points (xo, ^o) are the steady- 
state solutions to Eqs. and the ground state corre- 
sponds to the ones that give rise to the smallest energy. 
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FIG. 1: (Color online) Ground state molecular population j/o and 
energy gap A_E as functions of A. For (a) and (b), d = 0, the 
thin black lines represent the semiclassical results, dashed green 
lines are quantum results for N = 10, and thick red lines are 
quantum result for N = 100. For (c) and (d), d = 0.2 and only 
the semiclassical results are shown. 



Obviously x G [d, 1]. For convenience, we also introduce 
a variable 

y = 1 -x, 

which lies in the range of [0, 1 — d] and has the physical 
meaning that y/2 represents the normalized molecular 
population. For clarity, we will separately discuss the 
two cases: A = and A ^ 0. 



A. Case 1: A = 

In order to illustrate the effect of atomic population 
imbalance, we first present the results for d = 0. The 
ground state in this case is given by 



j/o = 1 , p undefined , 
y = i(VA J +3-A) 2 



</>o 



for A < -1 
for A > -1 



from which one can see that although yo is continuous 
throughout the A-space, the derivative dyo/dA has a dis- 
continuous jump at A = — 1. Therefore A = — 1 repre- 
sents a critical point that separates the pure molecule 
phase (yo = 1) from the atom- molecule mixture phase in 
the semiclassical theory. 

To study the corresponding quantum behavior and its 
connection with the semiclassical approach, we expand 
the Hamiltonian © using Fock state basis for a given 
set of N and D and diagonalize the resulting Hamilto- 
nian matrix. Both the quantum and the semiclassical 
results of ground state molecular population are shown 
in Fig. [lja). The quantum calculation always results in 
a smooth yo curve although it also shows a rapid change 
from to 1 in a small region near A = — 1. As expected, 
the quantum results approach the semiclassical limit as 
N increases. 

Further insights into the properties of the system can 
be gained by studying the excitations above the ground 



state. The quantum many-body excited states are ob- 
tained in the same manner as above through the diag- 
onalization of the Hamiltonian matrix. We are particu- 
larly interested in the "energy gap" , AE, defined as the 
energy difference between the first excited state and the 
ground state, which is plotted in Fig.QJb) for several dif- 
ferent N. The energy gap shows a minimum, which is al- 
ways finite, at the value of A around which yo rapidly ap- 
proaches 1. The semiclassical energy gap can be obtained 
through the following linearization procedure: Substitut- 
ing x — x + Sx and ip = ip + Sp into Eqs. (0} where 
(xo, Po) are the steady-state solution and (5x, Sip) repre- 
sent the small fluctuations away from the steady state, 
keeping terms up to first order in fluctuations, we have 

—Sx — 
dr 



dr 



Sp = < -A - 



y (1 - x o) (xq - d 2 ) cos p Sp, 
(1 - 3a;o) 



V(l - *o) - d 2 ) 
(d 2 + 2xo 



COS Po 



3x 2 ) 



A[(l-x a )(x 2 -d 2 )] 



3/2 



cos Po } & x > (6) 



where, in anticipation of later studies, we have not made 
the assumption of A = 0. The oscillation frequency of Sx 
and Sp can be derived straightforwardly as 



(d 2 + 2x Q 



4(1- 



,T 



(4 



d?) 



+ 3a;o-l 



cos tp 



x ) (xl - d 2 )cospo ■ 



(7) 



For ground state in the case of A = 0, the semiclassical 
excitation frequency reduces to 



y/A 2 + 3xo - 1 , 

which is the semiclassical energy gap. In particular, for 
d — 0, we have 



to 



x/A^ - !, 

[A 2 + 2-(VF+3-A) 2 /3] 



1/2 



for A < -1 
for A > -1 



which is plotted in Fig. QJb). The semiclassical energy 
gap vanishes at the critical point A = — 1 with a discon- 
tinuous jump in its derivative. 

Figure [Ha) an d (b) clearly show that the quantum 
result approaches the semiclassical limit as — > oo and 
hence the much simpler semiclassical theory is reliable for 
large N. Furthermore, there is a critical point at A = — 1 
for d = in the semiclassical theory which is absent in 
the quantum calculations with finite N, indicating the 
fact that no true quantum phase transition can occur in 
a finite system. 

We now discuss the case with finite atomic popula- 
tion imbalance, i.e., d ^ 0. Although semiclassical solu- 
tions to ground state population and excitation can be 
obtained analytically in the same fashion as in the pre- 
vious case for d — 0, the expressions are generally too 
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messy to be instructive. We therefore simply display the 
results in Fig. QJc) and (d) . Again we find that the semi- 
classical calculation reproduces the quantum result (not 
shown in the figure) in the large N limit. One major 
difference between d ^ and d — is that in the former 
there is no quantum phase transition even in the semi- 
classical limit: both the population and the energy gap 
changes smoothly as A varies, and the energy gap never 
becomes zero. 

From Fig. [TJ we can also see that starting from a pure 
two species atomic condensate, we can coherently create 
molecular condensate using the method of rapid adia- 
batic passage, e.g., by tuning A from a large positive 
value to a large negative value. Near perfect atom-to- 
molecule conversion [38] is achieved when A is swept adi- 
abatically ^39] which is confirmed by our numerical calcu- 
lations. However, as we demonstrate next, such a smooth 
conversion of atoms into molecules by a slow sweeping of 
A cannot be taken for granted when A ^ 0. 



B. Case 2: A^O 

With a finite A, the algebra becomes much more com- 
plicated. We resort to numerical calculations in this case. 
Consider first the semiclassical situation. The left panel 
of Fig. [5] illustrates the properties of the system with 
d = and A = — 5. Fig. [2] (a) and (b) show the molecu- 
lar population and mean-field energy for the semiclassical 
steady states. In the region A e [—3.56, — 1], there exist 
three steady states with similar energies as shown in the 
figure [4(J. The mean-field energy exhibits a swallowtail 
loop structure. Similar structures have been observed in 
condensates moving in optical lattice potentials [4l[ and 
in two-component condensates |42| under certain condi- 
tions, and are associated with dynamical instability. 

In our system, by calculating the excitation frequency 
using Eqs. ^ and ([7]), we find that one of the three 
steady states, represented by the red dashed lines in 
Fig. [2] (a) and (b), possesses imaginary excitation fre- 
quency, a signature of dynamical instability. This unsta- 
ble state links the two stable ones, representing a classical 
example of bistability which has been intensely studied in 
the context of nonlinear optics and laser theory [43| . The 
existence of such a state is the key to the development of 
atom-molecule switch, the matter wave analog [44j of op- 
tical bistable switch, for controlling matter waves by mat- 
ter waves in a coherent and bistable fashion. Under such 
a bistable situation, no matter how slow we tune A, the 
system will not be able to follow the ground state — when 
we enter the dynamical unstable region, a discontinuous 
jump will necessarily occur and the atom-molecule con- 
version efficiency will suffer. This is confirmed in our 
numerical simulation as shown in Fig. [21c) where we lin- 
early sweep A from a large positive to a large negative 
value starting from a pure two species atomic condensate. 
In this example, only about 60% of the initial atoms will 
associate into molecules. 
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FIG. 2: (Color online) Left panel shows an example of dynamical 
instability with d — and A = — 5. (a) Molecular population in 
the low lying semiclassical steady states as functions of A; the state 
represented by the red dashed curve is dynamically unstable, (b) 
Corresponding dimensionless semiclassical mean-field energies as 
calculated using Eq. J5j ; (c) Molecular population as A is linearly 
swept; (d) The corresponding quantum many-body energy spec- 
trum for N — 20, only the lowest five eigenenergies are shown. As 
the classical Hamiltonian JSJ represents energy per pair of atoms, 
the quantum eigenenergy (in units of G) has been rescaled by a 
factor of (iV/2) -1 . Right panel is the same as the left except for 
A = 5 where there is no dynamical instability. 



It is instructive to examine the situation from the 
quantum many-body point of view. Fig. [DJd) shows the 
five lowest eigenenergies of the quantum Hamiltonian (J3]) 
for N = 20. The quantum mechanical energy spectrum 
exhibits a net of anticrossings enveloped by a swallowtail 
loop structure that will morph into the semiclassical en- 
ergy diagram as shown in Fig.^b). Similar semiclassical- 
quantum correspondence was observed in two-component 
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FIG. 3: Stability phase diagram in the A-A parameter space. The 
black region are dynamically unstable, (a) d = 0; (b) d = 0.2. 



condensates 1451 . |46[ and condensates in double-well po- 
tentials [13, l4cf ! 

In comparison, the right panel of Fig. [2] shows a situ- 
ation without dynamical instability. In this case, rapid 
adiabatic passage results in a near perfect atom-molecule 
conversion, and the system follows the ground state 
closely as A is tuned. 

Figure [2] shows that in order to create molecular con- 
densate with high efficiency using the rapid adiabatic 
passage method, it is of crucial importance to avoid the 
unstable regimes Q . Fig. [3] shows the stability phase 
diagram in the A-A parameter space. We find that dy- 
namical instability occurs in the region of A < —1 and 
A < — 1 and is quite sensitive to atomic population im- 
balance d: With the increase of d, the unstable region 
shrinks. Therefore tuning the population imbalance pro- 
vides us with a handle to control the dynamical stability 
of the system. 



IV. COHERENT ATOM-MOLECULE 
POPULATION OSCILLATIONS 

Coherent po pul ation oscillation has been predicted 
[3 EE J3, 111, LUl H HJ HI and experimentally mea- 
sured 7| in systems of homonuclear molecules coupled 
to atomic condensate. Besides proving the phase coher- 
ence between atoms and molecules, a measurement of 
the oscillation frequency can tell us many properties of 
the system such as the molecular binding energy, atom- 
molecule coupling strength, etc. We therefore want to 
study in this section the population oscillation dynamics 
in our system starting from a pure atomic cloud, focusing 
again on the effect of atomic population imbalance. 



In a dissipationless system, the total energy is con- 
served so that the Hamiltonian represents another con- 
stant of motion and the semiclassical problem becomes 
integrable. For an initial state with pure atoms, i.e., 
i = l, the energy constant according to Eq. ([5]) is 
E = -A + A/2. By inserting 



Ax- 



COS(/3 



A 2 
2 x 



E 



- x) (x 2 - d 2 ) ' 



(8) 



which is obtained from Eq. (J5]), into Eqs. ([3]), we can 
easily find that 



+ (1 - A' A) y 2 - AV/4] , (9) 

where A' = A — A and y = 1 — x as before. 

The solution to Eq. © can be expressed in terms of 
the elliptical functions and strongly depends on the roots 
of the cubic equations inside the square bracket. A dis- 
cussion of the solution for the model with homonuclear 
molecules [d = 0) can be found in Refs. [l6|,|20||. Here, in 
order to gain physical insight into the effect of the popula- 
tion imbalance on the oscillation dynamics, we will focus 
on the simpler case with A = 0. Under this condition, 
Eq. © reduces to 



dy 



--Ay. 

4 y 



whose solution, when expressed in terms of Jacobi's el- 
liptic function, has the form 



y 



y- sn 2 [Jy^r/2, y/yl/y+j 



(10) 



where 



V- = n 



2/4 



4 



d 2 - 1 



A 2s 



4y- 



(11) 



Equation (|10[) describes an undamped oscillation in 
which y changes from to the peak value y- with a 
period 



4F 



T = 



Vy+ 



(12) 



where F(ir/2,k) is the complete elliptic integral of the 
first kind. 

We plot the amplitude j/_ and period T of the molecu- 
lar population oscillation with respect to A for different 
d in Fig. [2Ja) and (b), respectively. The figure is sym- 
metric with respect to A — so we only present the case 
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FIG. 4: (a) and (b), molecular population oscillation amplitude 
and period, respectively. The three curves correspond to d = 0, 0.2 
and 0.5 in descending order, (c) Molecular population dynamics 
for A = A = d — 0. Solid line: semiclassical result; dashed 
and dotted lines: quantum result for N = 100 and N = 1000, 
respectively, (d) On resonance semiclassical oscillation period as a 
function of d. 

with A > 0. From Eq. (fTTj) . we find that for any given 
d, the oscillation reaches a maximum value of 

V- = 1- d, 

at resonance, i.e., A = 0. 

One peculiarity from the semiclassical calculation is 
that when d = 0, the oscillation period diverges at A = 0. 
In this case we have y+ = y_ = 1 and Eq. [)10p becomes 

y = tanh 2 (r/2) , 

which shows that atomic (molecular) population de- 
creases (increases) monotonically until all the atoms are 
converted to molecules. The quantum mechanical cal- 
culation, however, does show damped population oscilla- 
tions under the same condition, as illustrated in Fig.HJc). 
The difference between the semiclassical and the quan- 
tum results arises because the former does not take atom- 
molecule entanglement into account. The same behav- 
ior will also occur in homonuclear molecule association 
and has been studied in Ref. In heteronuclear 

molecule association with finite d, the period T as given 



by Eq. (|12[) never diverges. Using the asymptotic formula 
for F(tv/2, k), one can show that, on resonance, 



VTTd d '■ 

for small d. The resonant oscillation period as a function 
of d is shown in Fig. [D(d). 

The situation becomes much more complicated in the 
case of A ^ and in general no simple analytic formula 
for population oscillation can be found. The general fea- 
tures are nevertheless still preserved: semiclassical result 
shows undamped oscillation while quantum calculation 
yields damped oscillation, and the quantum result ap- 
proaches the semiclassical limit as N increases. 



V. CONCLUSION 

In conclusion, we have studied coherent association of a 
two species atomic condensate into heteronuclear molec- 
ular condensate using a three-mode model, emphasizing 
the effect of atomic population imbalance. In particu- 
lar, the population imbalance, together with detuning 
and collisional interaction strength, will significantly af- 
fect the excitation and stability properties as well as co- 
herent population oscillations of the system. We have 
also carefully analyzed the differences and connections 
between the semiclassical and the quantum many-body 
treatments. 
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